Heavy quark asymmetries with DELPHI 
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The measurements of the forward-backward asymmetry in Z ^ cc and Z ^ bb decays are 
among the most precise determinations of sin^ ^^^j^j.. In this paper the results obtained by the 
DELPHI experiment at LEP with three different analyses are reviewed together with the impact of 
the combined LEP result on the global Electroweak fit. 



I. INTRODUCTION 

As a consequence of the parity violating couplings of the Z boson to the fermions, in electron positron annihilation 
at the Z mass fermions are more likely produced in the forward direction, with respect the incoming electron, than in 
the backward. In the Electroweak Standard Model the asymmetry at the Z pole |^ is expressed as: 
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Ae and Af are functions of the ratio a^/ of vector Vf and axial a/ couplings of the Z boson to the fermions: 
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The ratio depends on the quantum numbers of the fermion, the charge Qf and the weak isospin I3J, and on the 
fundamental parameter of the theory the electroweak mixing angle sin^ 9\y: Xf = l — 4:\Qf\ sin^ 6\y. 
The dependence of the Born level Af on sin^ 9w for different fermion species is shown in fig 0: For sin^ 6w 
the asymmetry for qq final states is from 3 to 5 times larger than for leptonic final states. 
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FIG. 1. Dependence of Born level ^/ on sin 9w for leptons, up-type quarks and down- type quarks. 

As among all the qq final states the decays of the Z boson into cc and bb pairs are the most easily distinguished, 
the experimental effort was concentrated on the measurement of Apg and Apg. 

Most of the electroweak radiative corrections can be accounted for by introducing in the Born level equations an 
effective electroweak mixing angle sin 0lY,eff defined for each fermion family. On the contrary the sensitivity of the 
quark asymmetry to the final state couplings is heavily suppressed. This can be understood from the the dependence 
of the sensitivity on and Xf: 



^ Throughout this paper pseudo-observables defined at the Z pole are indicated with an apex while values measured at the 
Z peak {y^ = 91.26 GeV) without. 
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The final state c and the h quarks have a smaller charge and a higher value of compared with the initial state 
leptons. The sensitivity to the quark couplings is reduced of about a factor 50 for the h and 5 for the c asymmetry. 

The main motivation for a precise measurement of the asymmetry is that m//, the mass of the still undiscovered 
Higgs boson, enters in the definition of sit? O'^ ^. Table ^ shows that, even if the dependence is only logarithmic, the 
measurement of Apg sets currently the most stringent limits on rriH'- A 4 % accuracy in the determination corresponds 
to values of the mass of the Higgs boson in a range from 70 GeV/c^ to 1000 GeV/c^. On the contrary the ratios of 
partial widths of the Z boson into cc and bb pairs are, in practice, not sensitive to ran at all. The current precision 
on Apg is smaller than this value and the prediction on the Higgs boson mass is spoiled mainly by the uncertainty 
on aQED{Ml). 

TABLE I. Dependence of the electroweak pseudo observables in the heavy flavours sector on radiative correction (from pi] 
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II. EXPERIMENTAL TECHNIQUES. 

All the analyses of the asymmetry follow three steps: The tag of the flavour of the decay of the Z boson, the 
determination of the direction of the primary quark and the separation of quark and antiquark hemispheres. The 
analyses described differ on how the tag is performed: Exploiting the long flight distance of the heavy flavoured 
hadron, the presence of a high momentum lepton or finally a reconstructed D meson in the final state 0. This choice 
determines also how the assignment of quark-antiquark hemispheres is done: Using the charge flow in the event or 
the charge of the identified lepton or the charge state of the D meson. Common to the analyses is the determination 
of the direction of the initial qq pair done with the thrust axis of the event conventionally oriented to form an angle 
9t < 90° with the direction of the incoming electron beam. The observed asymmetry is finally extracted from a fit 
to the differential asymmetry ^+_^^^~ where N^{Nf) is the number of forward (backward) events in bin i of cosOt- 

A. Apg using jet charge technique. 

Methods exploiting the long fly distance of 6-hadrons {'yf3cTb ^ 2 mm) and the accurate resolution on secondary 
vertex reconstruction of silicon vertex detectors provide the most efficient way to select a sample of bb events. In order 
to maximize the efficiency versus purity performance the 5-tag algorithm used in the DELPHI experiment combines 
the informations from the reconstructed secondary vertex (its effective mass, the rapidity of the tracks associated to 
it and the fraction of the jet energy carried by them) with the lifetime information. 



^If not explicitly mentioned charge conjugate states are implicitly included. 
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FIG. 2. Distributions of the 6-tag variable in data and simulation for 1994 data. 



In the measurement of Apg the fraction Pi,{cut) of bb events in the sample after a cut in the &-tag variable is 



determined from the data themselves using the relation: 

F{cut) - Rc X Scicut) - (1 



Pb{cut) 



Rc - Rb) X Eudsicut) 



F{cut) 

where F{cut) is the fraction of the data after the cut, £uds{cut) and £c{cut) the efficiencies for light {q — u, d, s) and 
charm quark taken from the simulation and Rc and Rb the partial widths of the Z into cc and bb pairs from the 
Standard Model. Fig. ^ shows the distribution of the 6-tag variable in data and in simulation. The value of the tag 
chosen for the analysis corresponds to et = 75% and Pb = 92%. 

The separation of b from b quark relies on the hemisphere charge Qhem defined as the sum of the charges qi of the 
tracks in each hemisphere, as defined by the thrust axis T, weighted by the projection of their momenta pi along T 
itself, to some power k: 

^hem — , . , -», 

\P^■T\- 

This estimator is based on the fact that, due to the electric charge conservation, the particles produced in the 
hadronization of the primary quark retain some information of its charge. In the analysis done by DELPHI k = 0.8. 
This choice minimizes the total statistical and systematic error. 

In each event the total charge Qtot — Qf + Qb and the charge flow Qfb ~ Qf ^ Qb are measured. A part from 
reinteractions (Qtot) — 0, while the average charge flow is directly related to the asymmetries: 

{Qfb) =Y.W^^FB (1) 

Pq is the fraction of events of flavour q in the sample, rjq is an acceptance correction factor and Sq ~ {Qq — Qq) is 
the charge separation: It would be twice the charge of the quark if quarks would be observed directly. As for P^, the 
charge separation for b quark is measured directly in the data. The principle of the method is sketched in fig. |^: For 
a pure b sample the charge separation leads to an increase of the spread apB of the distribution oi Qpb compared 
to utot from the distribution of Qtot so that 5^ ~ — c^qt- fo'" sample composition Su,d,s,c are taken 
from the simulation, carefully tuned in order to reproduce the measured distributions of hadronic event shapes and 
charged particle inclusive quantities. 
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FIG. 3. Sketch of the principle of the (Qfb) and the 5f measurement for a single (down type) flavour f. 
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The observable {Qfb) in data and simulation as a function of increasing &-purity of the sample is shown in fig The 
observed difference is due to different values of both and the input asymmetry between data and simulation. 
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FIG. 4. (Qfb) as obtained from data and simulation for k = 0.8 and for 1994 data. 

To exploit the angular dependence of the asymmetry itself, Apg is extracted from a fit of equation |l| to (Qfb) 
distribution in 4 bins of cos6't, the acceptance being limited to 9t > 35° by the actual angular coverage of the vertex 
detector. 

The results for the Z peak asymmetry for 1992-95 sample is Apg = 0.0982 ± 0.0047 (stat). The main contributions 
to the systematic error are listed in the top part of table |l[ The largest contribution (±0.0011) comes from the 
hemisphere-hemisphere charge correlation which is due to charge conservation, the common thrust axis and the 
particle crossovers between hemispheres. This correlation is checked between data and simulation by means of the 
variable H defined such that differences in the description of the secondary interactions mostly cancel out: H = 
{QfQb) — {QtotY I ^- The systematic error accounts for 20 % discrepancy between data and simulation found in 
the distribution of H. 

Effects related to the description of the physical processes which could effect the &-tag performance (gluon splitting 
into cc and bb, K" and A content in light quark events, lifetimes, spectra and production fraction of D mesons) are 
smaller than ±0.0003 each. 



TABLE II. Main contributions, in units of 10 ^, to the systematic error to Apg and Apg for Z peak data. 













Jet Charge tag 






r:; ± 5% 

H charge correlation 
Detector resolution 


± 0.4 
± 1.1 
± 0.7 




total 






± 1.6 


leptons tag 


c decay model 


T 1.8 




± 1.4 






± 1.9 


Px and jet reconstruction 


± 1.3 




bgd asymmetry 


± 4.7 


B°BO mixing 


± 1.1 




total 


± 6.3 




± 2.7 


D mesons tag 


MC stat 


± 2.5 




± 3.5 




fit method 


± 1.7 




=F 2.8 




bgd asymmetry (6, c quark) 


T 1.3 




T 3.6 

± 5.8 




total 


± 3.5 




=F 8.5 



B. a;:,'" and A^;^ using prompt leptons. 

Electrons and muons are produced in the decay of heavy flavoured hadrons mainly by 3 processes with a branching 
ratio of about 10 % each: 

• primary semileptonic h decays, b £~ 

• weak cascades of b hadrons, b c ^ £^ 

• primary semileptonic c decays, c ^ i~ 
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These leptons can be used to select on statistical basis cc or bb events as they have momentum (p) spectrum harder 
than fragmentation particles, {xe)c ~ 0.5 and {xE)b ~ 0.7, and high transverse momentum with respect the 
direction of the jet ^, because of the high mass of the parent hadron. 

The correlation between the charge of the lepton and the one of the quark is exploited to determine the quark 
hemisphere. The correlation is different in the three classes, resulting in a dilution of the measured asymmetry: 

kfi = {l-2x) {h - hc)A% - /cA^B + /bgd A^l'* (2) 

where the /i's are the composition of the sample, (1 — 2%) is the dilution factor due to the B^^B" mixing and the last 
term accounts for a non zero asymmetry of the misidentified hadrons (mainly punch through hadrons and photon 
conversions) and non prompt leptons from K and tt decays. 

DELPHI has recently introduced in the analysis of the data collected in 1994-95 the use of lifetime based 6-tag 
variable which is used to remove the decays of the Z into light quarks final states. This allows an estimation of the 
background level less dependent on the momentum of the lepton compared to the previous procedure. This consisted 
first in the selection of a high {p,p±) region, highly pure in the lepton content, which fixed the lepton identification 
efficiency and then in the evaluation of the background level in a low (jj,p±) region. Moreover the use of the 6-tag 
enriches the sample in cc events increasing the sensitivity of the measurement to Apg . On the contrary the reduction 
of the statistical error for the b asymmetry is smaller. This is because the useful sample is limited to the region of 
Pt > ^ GeV/c as the lower px region is equally populated of primary b and cascade leptons giving opposite sign 
contributions to the observed asymmetry: The 6-tag variable has little effect in separating these two classes of events. 

The rcsuhs for the Z peak asymmetry for 1991-95 statistics are Af-g = 0.0770 ± 0.0113 (stat) and A^^ = 0.0998 ± 
0.0065 (stat). The main contributions to the systematic error are listed in the central part of table ||. The various 
contributions can be separated into two categories: The terms arising from the model actually used to simulate the 
semileptonic decay processes and the ones related to the description of detector effects. In case of the b asymmetry 
they both amount at about ±0.002 the largest effect coming from the description of the c decay (±0.0014). 



C. Apg and Apg using reconstructed D mesons. 



A reconstructed D meson uniquely tags a hadronic decay of the Z boson into a heavy flavour qq pair as it can be 
produced only in the hadronization of a primarily produced cc pair or in the decay of a &-hadron in a bb event. The 
charge state of the reconstructed D meson is correlated with the charge of the parent quark. Therefore the values of 
Apg and Apg are extracted form a fit to the distribution of cos 9t in D meson events. 
The analysis done by DELPHI is based on the investigation of 8 decays channels: 

D*+ D07r+ 

- (K-7r+)7r+ 

— i- (K TT+TT 7r+)7r^ 

^ (K-7r+(vr«))^+ 
^ (K-M+z7^)^+ 
^ (K-e+;7,)^+ 

D" ^ K-7r+ 
DO ^K-7r+(7rO) 
D+ ^ (K-7r+)7r+ 

The first step in the analysis is the reconstruction of D" and D+ mesons. In case of D" a low momentum pion Tr^i 
[p* = 40MeV/c) with the correct charge correlation with the meson is searched for the D*+ reconstruction. The most 
relevant characteristics of the DELPHI detector for this part of the analysis are the particle identification provided 
by the RICH and the specific energy loss dE/dx measured by the TPC These informations are combined into a 



^The axis of the jet is defined excluding the lepton momentum. 

5 



pion veto used in the D*^/"*" channels to reduce the combinatorial background which mainly consists of misidentified 
pions. 

Further reduction of the background is achieved with cuts in the helicity angle of the kaon candidate with respect 
the D"/+ flight direction and in the 3D decay length of the meson. These cuts are applied in function of the scaled 
energy of the meson Xe{D) — 2Ey)I ^fs to account for the energy spectrum of charged particles in hadronic Z decays, 
peaked at low momentum. Finally candidates are selected either in a Am = mn* — tod" region for D*+ channels, or 
in a mass interval for a D°/+ decay. The range of Am values goes from 160 MeV/c^ to 250 MeV/c^ accordingly to 
the invisible energy in the final state, while for D+Z" channel a mass region ±200 MeV/c^ around the nominal mass 
is selected. The reconstructed mass spectra for two channels are shown in fig. |^ 
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FIG. 5. The mass difference Am spectra for D*^ (K t^^^t^^i and D*+ (K 7r^(7i"°))7rj; decay modes for data and 
simulation. 

To extract Apg and Apg the contributions of D mesons from cc and bb events should be separated. This is achieved 
fitting together Apg and Apg in bins of cos^t, scaled momentum Xe{D) and &-tag variable Vev (fig- 0)- 
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FIG. 6. The distributions of Xe{D) and &-tag probability for the signal region for D* 
simulation. 



(K 7r''')7r^ decay for data and 



In fig. H one can notice that, in addition to the pure combinatorial background, in the signal region there are two 
other components: Partially reconstructed D*+ mesons and reflections from other decays mode. Because of the charge 
correlation with the primary quark, these two categories have to be treated like signal in the fit. 

Finally in case of bb events an effective mixing is introduced to take into account the B'^B'' mixing for B ^ D 
produced mesons and the so-called "upper vertex" production of the charm in 6 ^ DD decays. This was computed 
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from data collected at the T(4s) resonance by the CLEO and the ARGUS collaborations. The values obtained are 
Xeff = 0.222 ± 0.033 for D+ and D*+ modes and Xeff = 0.170 ± 0.030 for the 0° channels, both different from the 
average x = 0.1214 ± 0.0043 measured at LEP for &-hadrons. 

At the Z peak about 62 x 10'^ D decays were reconstructed and the corresponding value of the asymmetries are 
A^B = 0.0659 ± 0.0094 (stat) and A^^ = 0.0762 ± 0.0194 (stat). The main contributions to the systematic error are 
listed in the bottom part of table Besides the limited MC statistics the main contributions are related to the fit 
method itself and to the residual asymmetry of D mesons from pure combinatorial background in cc and bb events. 



III. LEP COMBINED RESULTS. 

The currently available determinations of Apg and Apg from LEP experiments are shown in figure 0. The precision 
on the measurement of Apg is currently 7 % and the average is dominated by the measurements performed with the 
D meson tag. On the contrary, it should be noted that the most precise single determination comes from the lepton 
analysis tag from OPAL. This indicates that some improvements in the measurement of Apg can still be achieved. 
The precision on the measurement of Apg is currently 2 % with equal weights from jet charge and lepton tag analyses. 
The impact of these measurements on the determination of sin^ 0^ ^ is 
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sin^ e\v^^ff = 0.23225 ± 0.00038 for Ap'g 



Apg provides together Alr from SLD the most precise value available at the moment. On the contrary the two 
measurements are more than 2 sigmas apart from each other. 
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FIG. 7. 



Results for ASj'S and AI?^ shown at the ICHEP 



Vancouver Conference. 
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